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Abstract 

A system with Bose-Einstein condensate is considered in the frame of the self- 
consistent mean-field approximation, which is conserving, gapless, and applicable 
for arbitrary interaction strengths and temperatm'es. The main attention is paid 
to the thorough analysis of the condensate and superfluid fractions in the whole 
range of the interaction strength, between zero and infinity, and for all temper- 
atures between zero and the critical point Tc- The normal and the anomalous 
averages are shown to be of the same order for almost all interactions and tem- 
peratures, except the close vicinity of Tc. But even in the vicinity of the critical 
temperature, the anomalous average cannot be neglected, since only in the pres- 
ence of the latter the phase transition at Tc becomes of second order, as it should 
be. Increasing temperature influences the condensate and superfluid fractions in 
a similar way, by diminishing them. But their behavior with respect to the in- 
teraction strength is very different. For all temperatures, the superfluid fraction 
is larger than the condensate fraction. These coincide only at Tc or under zero 
interactions. For asymptotically strong interactions, the condensate is almost 
completely depleted, even at low temperatures, while the superfluid fraction can 
be close to one. 
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1 Introduction 



The relation between Bose-Einstein condensation (BEC) and superfluidity is of long- 
standing interest. The most thorough studies, both experimental and theoretical, of 
Bose-system properties have been accomplished for the region of low temperatures and 
weak interactions (see review works [1-7]), where the Bogolubov approximation [8,9] 
is applicable. In this region, almost the entire system is Bose-condensed, being just 
slightly depleted by interactions, at the same time, practically the whole system is 
superfiuid. When the condensate and superfluid fractions are so close to each other, 
one often uses as synonyms the terms "Bose-condensed" and "superfiuid". 

The opposite situation occurs for strongly interacting liquids, such as superfiuid 
^He, where the superfiuid fraction at low temperatures almost reaches one, while the 
condensate fraction never exceeds the values of the order of 10%. The superfiuid 
properties of helium are among the best-measured in experimental physics, as can be 
inferred, e.g., from the books [10-13]. The BEC fraction in superfiuid hehum, has 
been measured by using the x-ray scattering [14] and deep-inelastic neutron scattering 
[15]. Theoretical investigation, because of strong interactions between helium atoms, 
is rather complicated and mainly is done numerically, for instance, by means of Monte 
Carlo techniques [16,17]. 

It would be important to understand the behavior of the condensate and superfiuid 
fractions of the same system in the whole region of varying interaction strength and 
temperature. Present-day Feshbach-resonance techniques do allow for the variation of 
the interaction strength in a very wide range [18,19]. It is the aim of the present paper 
to investigate the condensate and superfiuid fractions for all temperatures between 
zero and the critical temperature and for all interaction strengths between zero 
and infinity. Such an analysis can be accomplished in the frame of the self-consistent 
mean-field theory [20-25], which is conserving, gapless, satisfies all thermodynamic 
relations and conservation laws. It was shown that this theory yields good agreement 
with Monte Carlo simulations for weak as well as strong interactions [24] and can also 
be applied for Bose systems in random potentials with arbitrary strong strength of 
disorder [25]. 

Studying here the condensate and superfiuid fractions, we analyse their properties 
both analytically and numerically. We use the system of units with fi = l and fc^ = 1- 

2 Uniform Equilibrium System 

Let us consider a uniform equilibrium Bose system with BEC. The appearance of 
BEC, as is known [26-28], is equivalent to the gauge symmetry breaking. The most 
convenient way to realize the latter is by means of the Bogolubov operator shift [29,30] 
representing the Bose field operator as the sum 



The first term here is the condensate wave function normalized to the number of 
condensed atoms 



'0(r) = ■q{Y) + iI)i{y) . 



(1) 




(2) 
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and the second term is the field operator of uncondensed atoms, satisfying the Bose 
commutation relations and the normalization to the number of uncondensed atoms 

Ni^J< A^(r)^i(r) > dr . (3) 

The angle brackets mean, as usual, statistical averaging. The condensate function 77 (r) 
and the operator of uncondensed atoms V'i(r) are treated as independent variables, 
orthogonal to each other, 

j 7;*(r)V^i(r) dr = 0. (4) 
The total average density is the sum 

P = Po + Pi (5) 
of the condensate density po and the density pi of uncondensed atoms, 

Po = — , pi = — , (6) 

where V is the system volume. The atoms are assumed to interact with each other 
through the local interaction potential 

$(r) = $o5(r) , $0 = 47r ^ , (7) 

m 

with Qs being the scattering length and m being atomic mass. 

For a uniform equilibrium system, the condensate function is constant, 

^(r) = ^/p~o. (8) 
The field operator of uncondensed atoms can be expanded over plane waves 



gik-r 

which gives 

^i(r) = • (9) 

k 

Prom the orthogonality equation (4) it follows that 

limafc = 0. (10) 

Therefore, in expansion (9), the summation is over k 7^ 0. This condition can either 
be shown explicitly in Eq. (9) or just one can keep in mind property (10). 

Accomplishing the Bogolubov shift for the field operator (1), we have the grand 
Hamiltonian 



ii-=^ii-W (11) 



n=0 
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consisting of five terms, labelled according to their order with respect to and a|.. 
The zero-order term does not contain the field operators of uncondensed atoms 

i/W = Q Po$o - /^o) iVo • (12) 
The first-order term is identically zero, 

H^^^ = , (13) 
because of the orthogonality equation (4) . The second-order term is 



k 

Respectively, one has the third-order term 

V 

P,Q 



+ 2po*o - l^ij alak + ^ Po*o (4«-fe + «-ik«fe) ■ (14) 



H^^^ ^ V t7 *o J2 («S«<z-p«p + alal_paq) (15) 
and the fourth-order term 

^^"^^ ^^12 aWqak+pag-k ■ (16) 

k,p,g 

In Eqs. (14), (15), and (16), the sums do not contain the terms with the operators a^, 
for which k = 0, due to the hmiting condition (10). 

The momentum distribution of uncondensed atoms is given by the normal average 

Uk = < alak > . (17) 

Because of the broken gauge symmetry, there also appears the anomalous average 

CTk = < akd-k > ■ (18) 

The normal and anomalous averages are equally important and neither of them can be 
neglected. Omitting the anomalous average would make the theory not self-consistent 
and would yield a spurious system instability [5,31]. Summing Eqs. (17) and (18) gives 
the density of uncondensed atoms 

Pi = ^ IZ^fc (19) 

k 

and the anomalous average 

<^i = 7E^fe' (20) 

k 



for which the value \ai\ defines the density of pair-correlated atoms 

Applying for terms (15) and (16) the Hartree-Fock-Bogolubov approximation and 
involving the Bogolubov canonical transformation 



ak = Ukhk + u*_f,blk , 



4 



we reduce [24] the grand Hamiltonian (11) to the form 



k 



in which 



is a nonoperator quantity; 



\ 



:2ii 



(22) 



(23) 



is the Bogolubov-type spectrum, but with the sound velocity c defined by the equation 



mc^ = (po + cri)^o ■ 
For the momentum distribution (17), we get 



"i^k^ 7— coth — - - , 



^k 



2T 



where T is temperature and 



ci^fc = 77— + mc^ 
2m 



And for the anomalous average (18), we find 

-,2 



We may notice that 



(7fe = - 



mc 
2^ 

<^k 



^k \ 



coth 

2Tj 



(24) 



(25) 



(26) 



(27) 



mc^ 2 

which indicates that au is, generally, of the same order as n^. 

Using Eq. (25) and passing in the standard way from summation over momenta to 
their integration, we obtain for the density of uncondensed atoms (19) 



Pi = 



{mcy 
"3^ 



1/2 



coth 



mc 

It 



X 



dx 



(28) 



For the anomalous average (20), using Eq. (27), we find 

1/2 



(7i = ctq - 



(mc) 



y/TTx^-l) 



2A/2 7r2 Jo 
where in the calculation of the term 

(Jo = 



coth 



mc 

It 



x \ -1 



(mc) 



dx , (29) 



(30) 



the dimensional regularization [24] is employed, in line with the general rules of the 
dimensional regularization [5,33]. 
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3 Condensate and Superfluid Fractions 

Our main aim here is to study the condensate and superfluid fractions. The condensate 
fraction 

no = — = 1 - ni (31) 
can be found by calculating the normal fraction 

-1 = ^ = -, (32) 

with the normal density (28). 

The superfluid fraction can be represented [2,25] as 



2Q 

ris = i - 

where Q is the dissipated heat, 



= 1 - , (33) 



expressed through the dispersion 

A2(P) = < P2 > - < P >2 

of the total momentum operator P. For the dissipated heat, in the considered mean- 
field approximation, we have 

Substituting here Eqs. (25) and (27), we get 

^^8^ I smh\ek/2T) (2^ ' ^^^^ 
The latter equation can be transformed to 

_ (mc)^ /"^ (a/T+x^ - 1)3/2 ^ 

A/2(27r) y/l + x^ sinh^(mc2x/2T) ' 

For the following analysis, it is convenient to introduce dimensionless quantities. 
We define the gas parameter 

7 = , (38) 

measuring the interaction strength, and the dimensionless temperature 

mT , , 

We introduce the dimensionless sound velocity 



and the dimensionless anomalous average 



(41) 



In terms of these notations, the dimensionless velocity (40), in view of Eq. (24), 
satisfies the equation 

= 47r7(no + a) . (42) 

The condensate fraction (31) is expressed through the normal fraction (32), for which 
we have 



Tlx 



37r2 



1/2 



coth 



'2t 



- 1 



dx 



The anomalous average (41), according to Eqs. (29) and (30), becomes 

.3 ,oo (v/rT^-i)V2 



2s2 



cr 



^3/2 2V2 
And the superfluid fraction (33) takes the form 

'■^ ( + a;2 - l^/'^x dx 



' re 
TT^ Jo 



coth 



' 2 ^ 
S X 



2t 



dx 



n. = 1 



gV^ttH Jo \^TTx^ smh'^{s^x/2t) ' 



(43) 



(44) 



(45) 



Equations (42) to (45), together with the relation uq — 1 — rii, define all characteristics 
we wish to investigate. 



4 Varying Interactions and Temperature 

Our aim is to study Eqs. (42) to (45) for the varying interaction strength 7 > and 
temperature t > 0. First, we find analytic expressions for low temperatures and for 
the temperature close to the critical point. 

A. Low Temperature 

At sufficiently low temperature, such that 

4 « 1 ' (46) 
we find from Eqs. (43) to (45) the asymptotic expansions for the normal fraction 

the anomalous average 

2s2 

V^o - TTT- , (48) 



TT' 



3/2 V / u -^2s 
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and for the superfluid fraction 



We may notice that between the condensate and superfluid fractions there is the rela- 
tion 

(1 -nj(l -no)^/^ ~ z — TTTT . (50) 

^ ^' 1357r(97r)V3 ^ ' 

Substituting these expansions into Eq. (42), we obtain the asymptotic behavior of the 
dimensionless sound velocity 

s ~ So + c£ {t^Q) . (51) 
Here the zero-temperature term is defined by the equation 

TT (475^ + Stts^ - 127r27)^ = 1927^5^ (Stt^ - sfj . (52) 
The coefficient a in the second term of Eq. (51) is 

a = 37r^(475g - Trgg - 127r^) 

16(27sg + 37rsg - 1877275^ + 3677^7) ' ^ ' 

Using expansion (51) in Eq. (48) gives the anomalous average 

(7 ~ (70 + (54) 

in which 



47sg + 37rs2 - I27r27 

^0 = TTT-n , = 



l + M(2853 + 7r2- 487r2a) 



• (55) 



127r27 ' 12^ 

Finally, we obtain the asymptotic temperature expansions for the condensate fraction 

and for the superfluid fraction 

~ 1 - (57) 
454 ^ ' 

Equations (56) and (57) show that Ug > uq. 

In order to specify the behavior of no and Ug as functions of the interaction strength, 
let us consider two limiting cases, of weak and strong interactions. When the interaction 
is weak, such that 7 — >• 0, Eqs. (52) and (53) give 

16 1 

So~20F7^/2^y 72, (7^0). (58) 
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Then the sound velocity (51) is 



s ~ 2 Vtt 7^/^ - — . (59) 

Remembering condition (46), we see that the considered expansions are vahd for the 
low temperatures for which 

t < 7 < 1 ■ (60) 
For the anomalous average (54), we get 

The condensate fraction (56) becomes 

g 3/2 / ^2 \ 

which is in agreement with the known temperature expansion for the weakly interacting 
Bose gas [5], first derived by Lee and Yang [34]. And for the super fluid fraction (57), 
we have ^ 

~ 1 - ( -] . (63) 



Noting that 

Us 873/2 / t2 

— ~ IH — IH 

riQ 30r \ 647 

we clearly see that > no, though they are close to each other when 7 <^ 1 and t <^1. 

In the opposite case of very strong interactions, when 7 — > 00, Eqs. (52) and (53) 
yield 

So ^ (377^) - ^ ( 3 j 7"^ " - ~ ^ (7 ^ 00) . (65) 
The sound velocity (51) behaves as 



For the anomalous average (54), we find 



4 ^7 Qtt 
In this way, we obtain the condensate fraction 



TT 



and the superfluid fraction 



no~^7-' + 0(t^) (68) 



2 / 57r \ 
^ ~ 1357r(97r)V3 + J " ^^^^ 



The ratio 



^ ~ — 7^ (7 ^ oo) (70) 



no TT 

shows that is much larger than Uq. 
B. Critical Region 

The analysis of Eqs. (42) to (45) demonstrates that there is the critical temperature 

where no, n^, a, and s all tend to zero. Temperature (71) is the same as the conden- 
sation temperature of the ideal Bose gas, as it should be in the case of a mean-field 
approximation for atoms with local interactions [5] . 
Considering the critical region close to tc, when 

f « 1 , (72) 
we find from Eqs. (43), (44), and (45) the normal fraction 

/ t \ s3 

"■-(^) (^^> 

the anomalous average 

2s2 St 

^^^V7^-2^' (74) 

and the superfluid fraction 

Calculating the last term in Eq. (75), the dimensional regularization was employed. 

When temperature t tends to tc, then it is convenient to introduce the relative 
temperature 

T = 1 - f ^ +0 , (76) 
which tends to zero. Then, solving Eq. (42) results in 



s ~ -y67r7 T^/^ (77) 
for any 7 > 0. Equations (73) and (74) yield the expansions 



ni ~ 1 - ^ T 2 W - 73 r^/^ 

2 V TT 



1/2 , , ,/22,^3/2 





tcr''^ + -^T+\l^tcr' 
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Using these, we obtain the condensate fraction 



(78) 



and the superfluid fraction 



Us 



As wc see, because of the ratio 



1 + 47r 



C(l/2)7 



amtc 



(79) 



no 



C(3/2)te 



(80) 



the superfluid fraction is again larger than the condensate fraction for all 7 > 0. 

Strictly speaking, the above expansions for the critical region are valid for not too 
large gas parameter 7. This is because the first term in Eq. (29) has been obtained 
using the dimensional regularization, which results in Eq. (30). The dimensional rcg- 
ularization is known [5] to be asymptotically exact in the limit of weak interactions. 
The analytical continuation of Eq. (30) to finite interactions may become not accu- 
rate for large values of the latter. In the above expansions, we constantly meet the 
combination 'jt, which is to be smaller than one to make the expansions quantitatively 
correct. From the inequality 'ft < 1, in the critical region, when t ~ tc, we get 7 < 0.3. 
Considering in what follows the critical properties of the system, for 7 > 0.4, we keep 
in mind that the behaviour of thermodynamic characteristics for these values of 7 is 
only approximate. 

The asymptotic behavior of Uq and Ug in expansions (78) and (79) makes it apparent 
that both riQ and tend simultaneously to zero at tc, so that BEC coincides with 
the superfluid transition. This phase transition is of second order, which agrees with 
the universality theory. According to the latter, the considered Bose system with 
BEC belongs to the universality class of the 3-dimensional 0(2)-symmetric spin model, 
hence, must display the second-order phase transition [5,35,36]. It is worth noting the 
importance of the anomalous average a. Though it tends to zero, as t — > tc, but it 
cannot be neglected, since no also tends to zero. Neglecting the anomalous average a 
would result in the first-order phase transition [5], which is not correct. In addition, 
omitting a would make the system unstable at all temperatures [31]. 

C. Numerical Solution 

In order to investigate the behavior of the characteristic quantities in the whole 
range of temperatures t e [0,tc] and for arbitrary interactions 7 > 0, we solve nu- 
merically the system of equations (42), (43), (44), and (45), together with the relation 
no — 1 — ni. Figure 1 presents the normal fraction (43) as a function of the dimen- 
sionless temperature (39) and the interaction strength (38). In Fig. 2, we show the 
anomalous average (44) as a function of the same variables t and 7. The condensate 
fraction is depicted in Fig. 3 and the superfluid fraction Fig. 4. Increasing temperature 
depletes both no and n^. But their dependence on the interaction strength 7 is not the 
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same. Increasing 7 always strengthens superfluidity, and Ug becomes larger. However 
the dependence of Uq on 7 is not trivial. At zero temperature, the condensate fraction 
diminishes with increasing 7. But at finite temperatures, no, first, increases with ris- 
ing 7 and then diminishes. This nonmonotonic behavior of no is connected with the 
nonmonotonic dependence of the anomalous average a on 7. The dimensionless sound 
velocity (40), given by the solution of Eq. (42), is displayed in Fig. 5. Temperature 
always diminishes s, while interactions make s larger. 



5 Discussion 

We have presented a detailed analysis of the condensate, no, and superfiuid, n^, frac- 
tions for arbitrary interactions strengths and for all temperatures in the internal [0, T^,]. 
The consideration is based on the self-consistent mean-field theory [20-25] . For the hm- 
iting cases of low temperatures and for those in the critical region, we derive analytic 
expressions for no and n^. And for the whole range of temperatures and interactions, 
we accomplish numerical calculations. The appearance of BEC and superfluidity oc- 
curs simultaneously at the critical temperature Tc. This transition is of second order, 
as it must be according to the universality theory. The superfiuid fraction is practically 
always larger than the condensate one. They coincide only at Tc or for the case of the 
ideal Bose case. 

It is important to stress the crucial role of the anomalous average (44). If it would be 
neglected, the phase transition would be of first order, which is not correct. And, more- 
over, neglecting the anomalous average renders the consideration not self-consistent and 
the system unstable. 

We may mention that the relation between the condensate and superfiuid frac- 
tions can be connected with the infrared behavior of the single-particle Green function 
Gii(k, ti;) in the momentum-energy representation. This asymptotic behavior, having 
the form 

|Gn(k,0)|:.^ (k^O), 

has been, first, obtained by Bogolubov, analyzed in great detail in Refs. [37], and 
summarized in his books [29,30]. As is evident, the Bogolubov infrared expression for 
G'ii(k, 0) can be rewritten as the hmit 

^ = lim^|Gn(k,0)| . 

The same asymptotic infrared expression has also been rederived and discused in Refs. 
[38-40]. Actually, the above relation is not an explicit definition of the ratio no/ug 
but it is an implicit equation, since Gn itself is a complicated function of no. In our 
paper, the relation between no and Ug is given by no = 1 — ni, with ni defined in Eq. 
(43), and by Eq. (45) for n^. These equations (43) and (45) also are not the explicit 
expressions for no and Ug, but are a part of the system of equations (42), (43), (44), 
and (45). Solving these equations makes it possible to extract the condensate and 
superfiuid fractions as functions of temperature and interaction strength. 

There is no simple general relation between the condensate and superfiuid fractions, 
represented as explicit functions of temperature and gas parameter. In some limiting 
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cases, it is possible to find asymptotic relations between these functions. For example, 
at low temperature, the fractions are related through Eq. (50). When, in addition, the 
interaction is asymptotically weak, satisfying Eq. (60), then both these fractions are 
close to each other, as in Eqs. (62), (63), and (64). 

In the opposite case, when temperatures are low, but the interaction is strong, the 
condensate fraction is drastically depleted, in agreement with Eq. (68). At the same 
time, the supcrfluid fraction, according to Eq. (69), is close to one. The limiting case, 
corresponding to Eq. (70), could be realized in low-tcmpcrature experiments with 
cold atoms by increasing their scattering length by means of the Feshbach-resonance 
techniques. Then the state of the Bose gas can be achieved, which contains a very 
small condensate fraction, though being almost completely superfiuid. 

Another way of creating the state with a tiny condensate fraction, but with a large 
superfiuid fraction, close to one, could be by loading atoms into an optical lattice, so 
that to reduce their effective mass. Diminishing the effective mass, as follows from Eq. 
(7), is equivalent to the increase of the effective interaction. In the one-dimensional 
case, the strengthening of atomic interactions would result in the effective fermioniza- 
tion of bosons, corresponding to the Girardeau mapping (see review article [7]). 

The analysis of the present paper can be straightforwardly extended to nonuniform 
systems, such as atomic systems in trapping potentials. For shallow traps, the local- 
density approximation (see Refs. [2.3]) can be employed. Then the overall consideration 
remains practically the same as above, merely with the appearing dependence on the 
real-space variable. In that case, the equations of the present paper can be interpreted 
as being written for the center of the trap. 

In the general case of an arbitrary trapping potential, we again could follow the same 
steps as above, just with some technical complications, caused by the system nonuni- 
formity [20-22, 32]. Then the main difference is that the Hartree-Fock-Bogolubov de- 
coupling for the field operators of uncondensed atoms should be done in the real-space 
representation and the general form of the Bogolubov canonical transformations [29, 30] 
has to be used, as in Ref. [32]. Following this way for atoms in a trap with a trapping 
potential U{r) and local interactions $(r) = $o(^(i"), we can use the Hartree-Fock- 
Bogolubov approximation for a nonuniform matter and the canonical transformations 
as in Ref. [32]. Then, introducing the notation 

^(r) = -^ + t/(r)+2$op(r)-/^i, 

A(r) = $o[po(r)+f7i(r)] , 

where the total density of atoms is the sum of the density of condensed and uncondensed 
atoms, 

p(r) =po(r)+Pi(r) , 

we obtain the Bogolubov equations, providing the diagonalization of the Hamiltonian, 

(jj{r)uk{r) + A(r)t;fc(r) = ekUk{r) , uj*{r)vkir) + A*{r)uk{r) = -£kVk{r) . 

Here k is a multi-index labelling the solutions to the above eigenproblcm. Since uj{r) 
is real, we can take A(r) to be also real. The Bogolubov equations define the spec- 
trum of collective excitations and the coefficient functions Uk{r) and Vk{r). The latter 
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functions have to satisfy the canonical conditions, similar to the uniform case, which 
would guarantee the Bose commutation relations of the field operators. Note that the 
anomalous average in these equations cannot be omitted, since, as has been shown 
above, the anomalous average is of the order of the normal average. Omitting the 
anomalous average would lead to incorrect results. 
The space-dependent condensate density 

Po(r) = \v{r)f 

is defined through the condensate wave function 77 (r), satisfying the equation 
V2 



2m 



+ U{r) 



?7(r) + $o + 2pi(r)^(r) + o-i(r)7y*(r) = Ato7y(r) 



If A(r) is real, then o'i(r) is also real. For an equilibrium system, the condensate 
function 77 (r) is real. The parameter /^o is given by the normalization condition 



^0 = / \vir)\' 



dr . 



The normal and anomalous averages, pi(r) and (Ji(r), are defined in the standard way. 
The number of uncondensed atoms is 



7Vi = j pi{v)dv 
The condition of the condensate existence 



min£fc = , £fc ^ , 

k 

gives 111 = fii(T) as a function of temperature and other system parameters, such as 
the gas parameter. Hence, the number of uncondensed atoms, A^i = Ni{T), is also a 
function of temperature and all other system parameters. As as result, we can find the 
number of condensed atoms A^o = No{T) as No{T) — N — iVi(T) and, respectively, the 
condensate fraction Nq/N. 

Thus, the whole procedure for trapped atoms is ideologically the same as for the 
uniform system. The technical complications come from the necessity to solve the 
equation for the condensate function and the Bogolubov equations for the coefficient 
functions and the collective spectrum. This, generally, requires the usage of extensive 
numerical calculations. 
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Figure Captions 



Fig. 1. Fraction of uncondensed atoms rii — ni{t,^) as a function of the dimen- 
sionless temperature t and of the interaction strength 7. 

Fig. 2. Anomalous average a = a{t, 7) as a function of the variables t and 7. 

Fig. 3. Condensate fraction uq — no{t, 7) as a function of the variables t and 7. 

Fig. 4. Superfluid fraction = ns{t, 7) as a function of the variables t and 7. 

Fig. 5. Dimensionless sound velocity s — s{t,'j) as a function of temperature t 
and the interaction strength 7. 
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Figure 1: Fraction of uncondensed atoms rii = ni{t,'y) as a function of the dimension- 
less temperature t and of the interaction strength 7. 
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Figure 2: Anomalous average a = cr(t,'j) as a function of the variables t and 7. 
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Figure 3: Condensate fraction no = no(t, 7) as a function of the variables t and 7. 
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Figure 5: Dimensionless sound velocity s = s{t,'j) as a function of temperature t and 
the interaction strength 7. 
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